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Given its weak coupling to bottom quarks and tau leptons, the Higgs boson 
may predominantly decay into invisible particles like gravitinos, neutralinos, or 
gravitons. We consider the manifestation of such an invisibly decaying Higgs bo- 
son in weak boson fusion at the CERN LHC. Distinctive kinematic distributions 
of the two quark jets of the signal as compared to Zjj and Wjj backgrounds 
allow to restrict the Higgs branching ratio to 'invisible' final states to some 13% 
with 10 fb~^ of data, provided events with two energetic forward jets of high 
dijet invariant mass and with substantial missing transverse momentum can be 
triggered efficiently. It is also possible to discover these particles with masses 
up to 480 GeV in weak boson fusion, at the 5a level, provided their invisible 
branching ratio is close to 100%. 



I. INTRODUCTION 

Some extensions of the Standard Model (SM) exhibit Higgs bosons which can decay into 
stable neutral weakly interacting particles, therefore giving rise to invisible final states. In 
supersymmetric models the Higgs bosons can decay with a large branching ratio into lightest 
neutralinos or gravitinos in some region of the parameter space [[Ip, leading to an invisible final 
state if R parity is conserved. Invisible Higgs decays also happen in models with an enlarged 
symmetry breaking sector, e.g. in Majoron models [0,0], where the Higgs disintegrates into 
light weakly interacting scalars. In the scenario of large extra dimensions, proposed by Arkani- 
Hamed, Dimopoulos and Dvali [Q], it is possible that the Higgs boson mixes with scalar fields 
arising from gravity propagating in extra dimensions [|^]. This mixing can lead to a sizeable 
invisible width for the Higgs. 

Certainly the presence of invisible decays modifies considerably the Higgs boson searches, 
making it much more difficult. At e+e^ colliders the problem is not very severe, and it has 
been shown that the Higgs parameter space can be probed completely [^. Presently, the LEP 
II collaborations exclude invisible Higgs masses up to 106.7 GeV [0]. On the other hand, 
the Higgs search at hadron colliders is much more difficult in the presence of such invisible 
decays. Previous studies have analyzed ZH and WH associated production ^ and ttH 
production [^ as promising channels. Considering statistical errors only and assuming that 
the Higgs couplings are the SM ones and that the invisible branching ratio effectively is 100%, 
the associated ZH production was estimated to be sensitive to Higgs masses niH ^150 GeV 
[P at the CERN LHC while tiH production might extend the Higgs mass range to 250 GeV 
[^. Both searches must deal with signal cross sections in the few fb range or below, however, 
and detector resolution and background normalization errors become crucial because signal 
and background distributions have virtually identical shapes after cuts. 

In this work we show that the LHC's potential to unravel the existence of invisibly decay- 
ing Higgs bosons can be considerably extended by studying Higgs production via weak boson 
fusion. Weak boson fusion (WBF) is a copious source of Higgs bosons at the LHC, yielding 
signal cross sections, after cuts, of order 100 fb, as we shall see. Furthermore, the presence 
of very energetic forward jets in the signal, with characteristic azimuthal angle correlations, 
allows us to tag the Higgs events and to efficiently suppress the backgrounds. It has been 
shown that WBF is among the most promising search channels in the 120-200 GeV Higgs 
mass range. It is possible to observe an intermediate mass Higgs boson through its decays 
^ ^ 77 &, H -> W*W* -^ e^fi^^T O, and even H -^ r+r' ^. Here we extend this 
list to include otherwise invisible decay modes. 

This paper is organized as follows. In Section II, we describe the techniques used to 
evaluate the relevant cross sections, while Section III contains the main characteristics of 
the signal and backgrounds and the cuts chosen to enhance the invisible Higgs signal. In 
Section IV, we present a detailed discussion of how to use data to predict the exact background 
normalization. Our results and conclusions are presented in Section V. 

II. CALCULATIONAL TOOLS 

We are considering the production of Higgs bosons in WBF, qq -^ qqVV -^ qqH {V = W 
or Z), with subsequent decay to undetectable particles. The signal is thus characterized by two 



quark jets, which typically enter in the forward and backward regions of the detector and are 
widely separated in pseudorapidity, and by a large transverse momentum imbalance, due to 
the Higgs' invisible decay products. Significant backgrounds can arise from any processes 
leading to two jets and missing transverse momentum, such as Zjj production with subse- 
quent decay Z ^> uu ox Wjj production with W ^ £u decay where the charged lepton is 
not identified. Another source are QCD dijet and multijet events with large missing transverse 
momentum generated by energy mismeasurements or from high px particles escaping through 
the beam-hole. These purely QCD processes possess a potentially very large cross section. In 
mismeasured two jet events, the missing px points in the direction of the jets and we use this 
fact to suppress this background. However, this is not true for the production of 3 or more jets 
and we include this QCD jjj background in our analysis. 

The signal and backgrounds are simulated at the parton level with full tree level matrix 
elements. This was accomplished by numerically evaluating helicity amplitudes for all sub- 
processes. Backgrounds include all order a^ real emission corrections to Drell-Yan production 
(to be called QCD Wjj and Zjj processes) and cross sections are calculated with code based 
on Ref. [|r3|]. The second large class of processes are the signal and the electroweak back- 



ground, EW Wjj and Zjj production, with contributions from the t-channel exchange of an 
electroweak boson off which the final state W, Z or Higgs is radiated. The code for these pro- 
cesses is based on Ref. [[T?|], and was checked with matrix elements generated using the pack- 



age Madgraph [jT^]. Madgraph code was also used to simulate the QCD jjj background at 
tree level. For all QCD effects, the running of the strong coupling constant is evaluated at one- 



loop order, with as{Mz) = 0.12. We employed CTEQ4L parton distribution functions [ [T^ ] 
throughout. The factorization scale was chosen as /i/ = mm(pT) of the defined jets. We took 
the electroweak parameters siii^ 6w = 0.23124, Oem = 1/128.93, mz = 91.189 GeV, and 
mw = 79.95 GeV, which was obtained imposing the tree level relation cos 9w = mw/mz. 
For the QCD jjj background, for which detector effects produce the missing px, we simulate 
experimental resolutions by smearing the energies (but not directions) of all final state partons 
with a Gaussian error given by A{E)/E = 0.5/v^ © 0.02 (E in GeV). 

W ^ lu decays in Wjj events lead to ajj;^T signature when the charged lepton I = e, ii,t 
is not identified. A precise modeling of this background requires a full detector simulation. 
We estimate the Wjj backgrounds by assuming that in the central region, \r]i\ < 2.5, all 
muons, electrons and taus with pr(0 > 5, 10, 20 GeV, respectively, can be vetoed, while any 
charged leptons below these thresholds will be misidentified and counted in the px balance 
only. In the forward regions, 2.5 < \r]i\ < 5, a lepton veto is taken to be impossible. Here, 
muons are assumed to give no pt deposit in the calorimeters, in contrast to electrons and taus 
whose entire energy is recorded. Note that the resulting Wjj background, within jet cuts 
given below, is about half the event rate of all Wjj, W ^ lu events with pr(z/) > 100 GeV, 
i.e. we are certain not to seriously underestimate the Wjj background. 

An important feature of the WBF signal is the absence of color exchange between the fi- 
nal state quarks, which leads to a depletion of gluon emission in the region between the two 
tagging jets. We can enhance the signal to background ratio by vetoing additional soft jet 



activity in the central region [|17[|. A central jet veto is ineffective against the EW Wjj and 
Zjj backgrounds which possess the same color structure as the signal. For the QCD back- 
grounds, however, there is color exchange in the t-channel and consequently a more abundant 
production of soft jets, with pt > 20 GeV, in the central region [jT^]. The probability of an 



event to survive such a central jet veto has been analyzed for various processes in Ref. [|Tq], 
from which we take the veto survival probabilities of Table | which are appropriate for the 
hard tagging jet cuts to be used below. 

The cross section for Higgs boson production via WBF is well known within the frame- 
work of the SM. We should keep in mind that this production cross section might be diluted 
in extensions of the SM. For instance, it is suppressed by factors sin^(/3 — a) or cos^(/5 — a) 
in supersymmetric models. Any suppression in the production cross section has the same ef- 
fect, for our study, as a branching ratio of invisible Higgs decays below unity, and we will not 
separate these effects in the following. 

III. SIGNAL AND BACKGROUND PROPERTIES 

The main features of the production of an invisible Higgs boson via WBF are the presence 
of two very energetic forward jets as well as a large missing transverse momentum. Therefore, 
we initially impose the following jet tagging cuts and missing momentum cut 

p^2.>40GeV , \Vj\<5-0 

\Vji - Vj2\ > 4.4 , r]ji ■ r]j2 < , (1) 

^T > 100 GeV . (2) 

A further reduction of the backgrounds, with good signal efficiency, is achieved by requiring 
a large invariant mass, Mjj, of the two tagging jets, 

Mjj > 1200 GeV , (3) 

and by selecting events where the azimuthal angle between the tagging jets, (pjj (measured in 
radians) is relatively small, 

0ii < 1 , (4) 

In order to motivate our choice of the |^y cut, we display, in Fig. [l], the ;^t spectrum after 
the cuts ([1]) and (Q), but without a central jet veto. The signal exhibits a peak around ;^t — 100 
GeV and it is much smaller than the backgrounds at small ;^t- Missing px generated by 
the QCD jjj background falls rapidly and this background becomes negligible above |^t = 
100 GeV. Note that we require (f)jj < 2.6 for the two tagging jets of the QCD jjj background, 
in order to avoid the soft singularities present near (pjj = n. Well above ;^t — 100 GeV, 
the missing pr spectra of the signal and the Zjj backgrounds have the same slope. Hence, a 
tightening of the ;^t cut soon becomes useless. 

The QCD backgrounds involve initial and final state gluons which tend to be softer than 
the quarks in WBF. In Fig. 0, this is reflected by the steeper falloff of the QCD backgrounds 
as Mjj, the dijet invariant mass, is increased. The Mjj > 1200 GeV requirement reduces 
these backgrounds sufficiently. Note that no central jet veto is included in Fig. |[ A further 
improvement of the signal to background ratio is possible by tightening the Mjj cut, but this 
will not be pursued in the following. 

The most distinct remaining difference between the Higgs signal and all backgrounds is 
the azimuthal angle correlation of the two tagging jets. The (f)jj distributions within the cuts 



(|T]-H) and including the central jet veto efficiencies of Table | are shown in Fig. ^j. The hV^V^ 
coupling of the Higgs boson in WBF favors Higgs emission opposite to both tagging jets, 
which leads to small values of (pjj. The Wjj and Zjj backgrounds, on the other hand, are 
smallest in this region and prefer back to back jets. The distinct shape difference of the (pjj 
distributions provides very powerful tools to probe for a Higgs contribution, either by a full 
shape analysis or by the (f)jj < 1 cut of @. 

The effect of such a cut is presented in Tables ^ and 0, where background and signal 
cross sections are given after imposing the cuts (^^ and (|l]-0). With an integrated luminosity 
of 10 fb^^ a total of 1670 background and of order 400 to 1000 signal events with (^^^ < 1 
are expected for an invisible branching fraction Br{H -^ invisible) = 1, giving a highly 
significant signal when statistical errors only are considered. 

rV. PREDICTING THE BACKGROUND 

Finding, or constraining, an invisibly decaying Higgs boson signal in jj;^T events is essen- 
tially a counting experiment since a resonance in the invariant mass distribution of the Higgs 
decay products cannot be extracted. The sensitivity of the search is thus determined by the 
precision with which the background rate in the search region can be predicted. Since the 
signal selection is demanding, including double forward jet tagging and central jet vetoing 
techniques whose acceptance cannot be calculated with sufficient precision in perturbative 
QCD, the background levels need to be determined directly from LHC data. 

Fortunately, a sizable sample of Vjj events {V = W or Z), with fully identified charged 
leptons from the V decay, will be available within the hadronic acceptance cuts discussed 
above. For Zjj events with two identified charged leptons, the |^t > 100 GeV cut is equivalent 
to Pt{Z) > 100 GeV and the only difference to the Zjj, Z ^ vv background is due to 
minimal charged lepton p^ and rapidity cuts which are needed to insure their observability. 
The cross section for QCD and EW Zjj, Z -^ l^l^ events (/ = e, /i), with pril) > 15 GeV, 
\r]i\ < 2.5 and within the cuts (|l]-^, including the central jet veto probabilities of Table |, is 
about 87 fb, of which 16.6 fb are expected in the (pjj < 1 region. Thus an effective luminosity 
of 10 fb^^ is sufficient to measure the Zjj background in the two regions (pjj > 1 and (f)jj < 1 
with relative errors of 3.8% and 7.8%, respectively. 

The level of the Wjj background can be extracted from W^ —^l^u events, with hadronic 
cuts as given above and the additional requirements priW) > 100 GeV and Pt(0 > 25 GeV, 
\r]i\ < 2.5. The combined cross section for QCD and EW {W -^ l^)jj events, including 
central jet veto, is 184 fb (975 fb) with (without) the (pjj < 1 requirement. With an effective 
luminosity of 10 fb^^ this translates into a statistical error of a mere 2.3% for the prediction 
of the total Wjj background within the cuts ([l]-0). Central jet veto efficiencies will be almost 
identical for charged lepton and missing px signatures of the decaying Ws and Zs. Other 
systematic errors like luminosity uncertainties, jet reconstruction efficiencies, or knowledge of 
the fractional contributions from WBF and QCD production are also eliminated by obtaining 
the background normalization from the leptonic W, Z data. Even trigger efficiencies can be 
determined directly from these events. 

The background normalization error can be further reduced by expanding the calibration 
region, at the price, however, of introducing a QCD uncertainty, due to the necessary extrap- 
olation to the signal region. For the Zjj background in particular, one would like to use the 



entire (pjj region to determine the background normalization, since the data sample within (Q) 
will be relatively small. These extrapolation uncertainties are small, as we will now show for 
the shape of the (f)jj distribution. 

At present we only have leading order (LO) calculations of the Wjj and Zjj QCD back- 
grounds available. Due to the small difference in weak boson mass, as compared to e.g. 
the large dijet mass required in our event selection, QCD corrections for these processes are 
expected to be very similar and we only analyze the shape of the (pjj distribution for the 
QCD Zjj background in the following. Shown in Fig. ^a) is da/dcpjj for four different 

renormalization scale choices, /i^j = y -5/4 (dash-dotted line) where s is the squared parton 



center-of-mass energy, /i^j = J{E^{Z) + Ptj^ +Ptj2)/^ (dashed line), the default choice 
'^s(/^?j) = '^siPTji)cisiPTj2) (solid line), and /x^j = Et{Z) (dotted curve). One finds that the 
normalization of the QCD background changes by up to a factor of 3 between these choices 
and another variation by a factor of 3 to 4 is obtained by changing individual renormalization 
scales between jir = yU/^/lO and fXR = lO/i^j- While the normalization of the QCD Zjj cross 
section changes drastically, the shape of the (pjj distribution is essentially unaffected. As a 
measure of shape changes we plot the fraction of events with (pjj < 1, 

as a function of ^, the scale factor for the four different renormalization scale choices jir = 
^/jPj^ listed above. The ^ dependence shown in Fig. ^b) is very small for individual choices of 
H% smaller in fact, than the differences between the four basic scales fi^. 

Fig. |(b) indicates that Ri = 0.19 ± 0.02 at LO QCD, i.e. a scale uncertainty of at most 
10% is found for the shape of the (pjj distribution, a very small QCD error indeed for a leading 
order calculation. While this error is still uncomfortably large for the actual determination of 
an invisible Higgs contribution to jj;^T events, a NLO calculation, which hopefully will be 
available by the time the experiment is performed, should push the QCD shape uncertainty to 
well below the 5% level. We assume a 5% QCD uncertainty on Ri in the following. 

Compared to this extrapolation from opposite side to same side dijet events, minimal 
changes in hadronic event properties are encountered when extrapolating from Vjj events 
with observed leptons, as discussed above, to the corresponding jj;^T sample. Hence QCD 
uncertainties for ratios like a{Zjj, Z — > uh')/a{Zjj, Z — > l^l^), within analog cuts, are ex- 
pected to be small. Indeed, an analysis of the scale uncertainties of these ratios, analogous to 
the one performed above for Ri, points to LO QCD errors of 3% or less for the Wjj ratios 
and 1-2% for the Zjj ratios. At NLO these QCD uncertainties should be entirely negligible 
and we do not consider them in the following. 



V. DISCUSSION 

The previous discussion points to at least two methods for constraining, or discovering, 
an invisible Higgs decay channel. The first uses the shape of the (pjj distribution only. With 
10 fb^^ already, the estimated 1670 background events in the (pjj < 1 region imply that Ri 
in (^ can be measured with a statistical error of 2.2%, which is already small compared to 



the systematic error of 5% due to QCD scale uncertainties. Taking the expectations for an 
Mh = 120 GeV Higgs boson as an example (see Table III), an invisible effective branching 
fraction Bh = Br(/f -^ invisible) would give rise to a modification 

1670 + 967a„^ 
9180 + 2380Bh ^ i 

In absence of a signal, a deviation of Ri of more than 10.7% can be ruled out at 95% CL, 
which translates into an upper bound 

Bh<M% at95%CL, (7) 

from shape information alone. Clearly this bound is systematics dominated and can be im- 
proved if it can be shown that our estimate for the QCD scale uncertainty is too conservative. 
A much higher sensitivity is achieved by making use of charged lepton signatures for 
Vjj events, as discussed in the previous section. Again assuming an integrated luminosity of 
10 fb^\ the Wjj background can directly be determined from W ^ lu events in the region 
(pjj < 1, with a statistical error of 2.3%. The number of Zjj events with Z -^ l^l~ , in 
this region, is still modest, which makes it worthwhile to combine the direct determination 
with an extrapolation of the {Z — > l^l~)jj cross section in the (pjj > 1 region, in spite of 
the 5% systematic uncertainty of Ri. Combining all errors in quadrature, the {Z -^ i^i^)jj 
background in (pjj < 1 can be predicted, directly from the data, with a 5.3% error. The 
combined background from Wjj and Zjj sources can thus be predicted with an accuracy of 
3.0%. Combined with the expected statistical uncertainty of the 1670 background events, one 
finds that the Higgs contribution can be measured with a la error of 6.4 fb, which, in absence 
of a signal, translates into a 95% CL bound of 

Ba{qq -^ qqH, H -* invisible) < 12.5 fb. (8) 

For an integrated luminosity of 100 fb~^ the combined backgrounds can be obtained with a 



precision of 1.2% and the bound (||) reduces to 4.8 fb. We present in Table |V| the invisible 
Higgs branching ratio that can be probed at 95% CL as a function of Mh- In the event that 
the invisible Higgs production cross section is suppressed with respect to the SM, the limits in 
this table apply to the product of the invisible branching ratio times the suppression factor. We 
find that the WBF channel extends considerably the sensitivity for invisibly decaying Higgs 
bosons, with respect to the ZH and itH channels [|J|]. With 10 fb^^ (100 fb^^) of data it 
is possible to discover these particles with masses up to 480 GeV (770 GeV), at the 5cr level, 
provided their invisible branching ratio is 1 . 

A disadvantage of the WBF process, as compared to ZH and itH production, is the re- 
quirement of a more complicated trigger. The latter provide hard isolated leptons as effective 
trigger signatures while the WBF signal possesses a dijet and missing px signature only. The 
two jets of the WBF signal are very hard and widely separated, however, which allows for the 
setup of a specific WBF trigger, with modest reliance on the Higgs decay signature. The QCD 
dijet cross section for events with two jets of pt > 40 (20) GeV within \r]j\ < 5, separated by 
\Vji ~ Vjil > 4.4, is about 2.3 (41) jJo. Requiring the presence of |^t with more than 3 to 4a 
significance, a condition which is easily fulfilled by the signal events with their ^2- > 100 GeV, 
will immediately lower these cross sections to the low nb level, i.e. to event rates of a few Hz 



only at a luminosity of 10^^cm~^sec~^. Similarly, the presence of isolated electrons, muons or 
photons, of modest pt, will allow for efficient triggers of H ^ tt, H -^ WW and iJ — i> 77 
events in WBF [p^|-|T^]. We strongly urge the LHC detector collaborations to develop such 



triggers for a broad-based study of WBF processes. 
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signal Hjj QCD Zjj and Wjj EW Zjj and Wjj 


Psurv 0.87 0.28 0.82 



TABLE I. Survival probabilities for the signal and background for a veto of central jets with 



PT > 20 GeV. From Ref. [18]. 





QCD Zjj 


QCD Wjj 


EW Zjj 


EW Wjj 


total 


a 


1254 


1284 


151 


101 


2790 


^surv 0" 


351 


360 


124 


83 


918 


Psurv (^{4'jj < 1) 


71.8 


70.2 


14.8 


9.9 


167 



TABLE II. Total cross sections (in fb) for the backgrounds after applying the cuts (|l]-|l) (first two 
lines) and <Q). In the last two lines we also include the central jet veto survival probabilities of Table |. 



Mh (GeV) 


110 


120 


130 


150 


200 


300 


400 


a 


282. 


274. 


266. 


251. 


214. 


154. 


110. 


^surv ^ 


245. 


238. 


232. 


218. 


186. 


134. 


95.7 


Psurv Cr{(l)jj < 1) 


99.4 


96.7 


94.3 


89.2 


77.0 


56.3 


40.7 



TABLE III. Same as Table || for the signal at several invisible Higgs masses, assuming 
Bt(H —>■ invisible) = 1. Cross sections are given in fb. 



Mh (GeV) 


110 


120 


130 


150 


200 


300 


400 


lOfb-i 


12.6% 


13.0% 


13.3% 


14.1% 


16.3% 


22.3% 


30.8% 


lOOfb-i 


4.8% 


4.9% 


5.1% 


5.3% 


6.2% 


8.5% 


11.7% 



TABLE IV. Invisible branching ratio that can be probed at 95% CL as a function of Mh for an 
integrated luminosity of 10 fb^^ and 100 fb^^. A SM production cross section is assumed. 
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FIG. 1. Missing transverse momentum spectra within the cuts (|l]) and (||). Results are shown 
separately for the EW Zjj (blue dashed line) and W jj (blue dotted line) backgrounds, as well as 
the QCD processes Zjj (black dashed line), Wjj (black dotted line), and jjj (magenta histogram) 
production. We exhibit the invisible Higgs contribution for Mu = 120 (red solid line) and 300 GeV 
(red dot-dashed line). 
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FIG. 2. Dijet invariant mass distributions when applying the cuts of Eqs. ^^. The lines follow 
the same convention as in Fig. |l|. 
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various background processes and the Higgs signal at Mh = 120 and 300 GeV. Results are shown 
after applying the cuts ([T}]3|) and including the effect of a central jet veto with the survival probabilities 
of Table |. The lines follow the same convention as in Fig. |l|. 
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FIG. 4. Scale dependence of the shape of the dijet azimuthal angle distribution of QCD Zjj events 
with Z ^ vv m LO QCD. Results are shown for four choices of the basic renormalization scale, 

//O - 



^/s/'i (dash-dotted Une), /i^j = \/{E^{Z) +Ptji ^Ptj2^/^ (green dashed line), the default 
choice a^(//jj) = as{pTji) "s(PTj2) (blue solid line), and yP^ = Et{Z) (red dotted curve). In (b) the 
fraction of events with (/>,, < 1 is shown as a function of ^, where /i/j = ^/Up- 
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